Acceleration of chemical reaction by chaotic mixing 
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Theory of fast binary chemical reaction, A + B — > C, in a statistically stationary chaotic flow 
at large Schmidt number Sc and large Damkohler number Da is developed. For stoichiometric 
condition we identify subsequent stages of the chemical reaction. The first stage corresponds to the 
exponential decay, oc exp(— At) (where A is the Lyapunov exponent of the flow), of the chemicals in 
the bulk part of the flow. The second and the third stages are related to the chemicals remaining 
in the boundary region. During the second stage the amounts of A and B decay oc 1/yi, whereas 
the decay law during the third stage is exponential, oc exp(— ■yt), where 7 ~ X/ySc. 

PACS numbers: 47.70.Fw, 47.27.Qb 



A common expectation is that random advection 
should essentially accelerate chemical reactions in fluid 
phase, since it should lead to homogenization of the re- 
action mixtures. Then dynamics is determined by an 
interplay of three factors: diffusion, advection and chem- 
ical reaction rate. Typical situation realizing in chemical 
reactors is that chemical reaction itself is much faster 
than mixing and diffusion, i.e. the Damkohler number, 
Da, which is defined as the ratio of the mixing time to 
the characteristic time of the reaction 0], is large. For 
the binary reaction this separation of temporal scales re- 
sults in formation of lamellar structure, build of stripes, 
populated solely by one chemical. The stripes of differ- 
ent chemicals are separated from each other by an in- 
terface of complicated shape, and the chemicals co-exist 
only in the narrow interface domain where the chemical 
reaction occurs. The reaction is limited by diffusion in 
the sense that diffusion controls fluxes of the chemicals 
into the interfacial reaction zone The physical 

picture of the acceleration due to the random advection 
is that it stretches domains populated by one chemical 
into thin sheets, so that the chemical reaction driven by 
diffusion proceeds more efficiently because of an essential 
increase of the interface area. In this letter we explain 
how this general physical picture, formulated initially for 
unbounded flows applies to chaotic flows confined to a 
finite geometry. 

We consider a binary chemical reaction, A + B — ► C, 
in a dilute solution of two chemicals. We study decay 
problem, with an initial distribution of the chemicals A 
and B, created by injecting solution of one chemical, say 
of A, into solution of the other chemical, B. It is assumed 
that the inverse reaction C — > A + B is negligible, i.e. 
there is no back influence of C on the distribution of A 
and B. Then molecular concentrations of the chemicals, 
n a and rif,, vary according to the following non- linear 
governing equations (see, e.g., 0) 

d t n a fi + (v ■ V)ra a ,f, = K a ^V 2 n a ^ b - Rn a n b , (1) 

where R is the reaction rate coefficient, v is the velocity of 



the flow, which is assumed to be incompressible (V • v = 
0), K a> b are the diffusion coefficients of the chemicals. 
Here, we assume that the fluid dynamics is independent 
of the chemical reaction, that is the velocity field does not 
sense changes in the chemical concentrations nor heat 
is released in the result of the chemical reaction. Our 
approach is also applicable to the case, realized in tubular 
chemical reactors, when the solution of the chemicals, 
prepared at the entrance, is then pushed through a pipe. 
In this case the position along the pipe plays the role of 
time in the decay problem. 

The major question addressed in the letter is: how do 
the total amounts of chemicals, N a ^ = J dr n a ^(t, r), 
decay as time t advances? We focus primarily on the 
stoichiometric case N a = Nf,. This case is of major in- 
terest for applications, as it allows to get pure product 
C (not mixed with the reagents), by the time reaction 
is completed. (An effect of a mismatch between N a and 
Nf, is briefly discussed later in this text.) We identify 
major stages of the chemical reaction and relate them to 
the chemicals decay in different parts of the flow. An 
essential part of the evolution is related to the boundary 
region. 

We discuss mainly the case of n a = Kf, = k. (It is 
argued later in the text that n a ^ Kb does not lead to 
significant changes in the theory.) Then one obtains a 
closed equation for the difference field, n = n a — rib, 



dtn + (v ■ V)rt = kV 2 ?i , 



(2) 



from Eq. (JIJ, i.e. one finds that n(t,r) is a passive 
scalar field. Note, that n has no definite sign, and that 
/ dr n = in the case of perfect matching, N a = A5. 

We assume that the chaotic statistically steady velocity 
field v(t, r) contains only few harmonics of the reservoir 
size, i.e. the flow is smooth. This regime can be realized 
in chemical reactors with mechanically rotating mixers 
or externally driven magnets stirring the fluid in the per- 
fect mixing devices and also in the tubular reactors at 
moderate Reynolds numbers. (See for a discussion of 
the chemical engineering principles behind various reac- 
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tor designs.) Complementary to its practical significance, 
passive scalar advection in a smooth chaotic flow is also 
a well studied (by both theoretical S H H E3 and 
experimental \v& Il3l Hif means) subfield of statistical 
hydrodynamics. (See also reviews FTgI The pas- 
sive scalar decay theory, developed in |llLll7| for an un- 
bounded flow, was recently modified for bounded flows, 
i.e. for chaotic flows with suitable (no slip) conditions 
on the boundary Smoothness of the flow allows 

one to approximate the velocity difference between close 
points by a linear, although fluctuating in time, profile. 
In the bulk region, the linear profile approximation is 
valid for separations smaller than the system size L. In 
the periphery, i.e. close to the solid boundary (wall), the 
linear profile approximation is valid for velocity fluctua- 
tions on a scale smaller than distance to the boundary. 
An important consequence of the linear velocity profile 
approximation is that close Lagrangian trajectories di- 
verge exponentially in time. The mean logarithmic rate 
of the nearby Lagrangian trajectories divergence defines 
the Lyapunov exponent of the flow, A. Notice, that in the 
peripheral domain advection is essentially anisotropic, 
and the stretching rate along the boundary is estimated 
by A, while the stretching rate in the direction normal to 
the boundary is significantly smaller. 

We now discuss characteristic spatial scales in the 
problem. The size of the system, L, which is also the 
chaotic flow typical eddy scale, is the largest scale in the 
problem. A comparison of the advection and diffusion 
terms in Eq. Q sets the dissipative scale of the flow, 
which in the bulk region is = \J n/\. We assume 
that the Schmidt number, Sc ~ (L/r^) 2 , is large, i.e. in 
the asymptotically wide range of scales, L ^S> r ^S> rd, 
advection dominates diffusion. The width of the diffu- 

1 /4 

sive boundary layer is estimated by ru ~ Sc ' rd, i.e. 
Tu > rd- Yet another important scale, associated with 
the chemical reaction itself, is the size of the reaction 
zone, r c h (the width of the interfacial domain where the 
chemical reaction occurs). In the bulk region the scale 
is estimated by r c h — rd[X/(Rn m )) 1 ' 3 , where Tiyji is a 
typical concentration of the chemicals inside the layers. 
(The estimate for the width of the reaction zone should 
be modified near the boundary, where it appears to be 
larger than in the bulk.) Initially, r c h is much smaller 
than rd (and, consequently, than L); the inequality is a 
consequence of the Da 3> 1 assumption. (Indeed, in ac- 
cordance with the definition, the Damkohler number can 
be estimated as Da ~ Rn /X, where no is a typical value 
of the initial chemical concentration. Thus, at t = 0, 
r c h ~ Da -1 / 3 ?^.) However, r c h grows as n m decreases. 
Thus even though the separation of scales is perfect ini- 
tially, it eventually breaks down at the latest stage of the 
chemical reaction. A cartoon illustration of the scale hi- 
erarchy is shown in Fig. ^ The magnified striped struc- 
ture is shown on the chart in the upper right corner of 
the figure. Regions populated by one chemical are single- 
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FIG. 1: Chemical reactor boundary is drawn by solid line on 
the main chart. Dotted lines separate bulk, peripheral and 
boundary domains of the flow. Chart in the upper right corner 
shows schematically the stripe structure magnified from the 
bulk and/or peripheral domains. Black and white regions are 
the ones populated by A and B respectively. The chart in the 
low right corner shows (under even stronger magnified glass) 
distribution of chemicals normal to their contact interface. 

colored. To resolve the interface domain, even stronger 
magnification is needed. Dependence of the chemicals 
concentrations on the coordinate normal to the interface 
is shown schematically on the chart in the lower right 
corner of the figure. 

The separation of scales, r c h <C rd, L, allows an im- 
portant simplification in the description. Indeed, the 
chemical reaction takes place in the r^-narrow interface 
domain, where the values of n a and rib are comparable. 
Outside this narrow region, i.e. in the region dominated 
by one of the chemicals, the presence of the other chem- 
ical is negligible. Thus, in the limit r c h — > 0, i.e. when 
the reaction zone becomes infinitesimally thin, one ob- 
tains n a — n, Tifc = for n > and n a = 0, rib = — n 
for n < 0. These relations imply a remarkable conclu- 
sion 0, Q : the fast chemical reaction can be described in 
terms of the linear setting (J2J) which does not contain the 
chemical reaction rate coefficient R. The reaction rate 
is determined by the diffusion fluxes of A and B to the 
n = surface. These fluxes are equal to each other and 
opposite in sign, which is translated, at r c h — > 0, into a 
continuity condition for Vn at the interface. This obser- 
vation also means that, while r c h is much smaller than 
all other relevant scales, our problem is reduced to the 
problem of scalar decay in chaotic bounded flow. 

Furthermore, from Eq. it is straightforward to de- 
rive equations for correlation functions of n (the deriva- 
tion procedure is similar to the one described in (IcUlU l. 
The equation for the mean value of n, (n), derived by av- 
eraging over times larger than the correlation time of the 
flow ro, is 

d t (n) = V Q (D a0 V/5 (n) ) + kV 2 (n) . (3) 
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In the case of a short-correlated (in time) flow, At <C 1, 
the turbulent diffusion tensor D is expressed through 
the velocity pair correlation function: D a p(r) = 
J °° dt(v a (t,r)v@(0, r)). Eq. remains valid in the 
general, not necessarily short-correlated limit. However, 
the relation between the eddy-diffusivity tensor and the 
velocity correlations becomes more complicated. D(r) 
tends to zero when r approaches the boundary, since the 
flow velocity tends to zero there. (The longitudinal com- 
ponent £)ii of the tensor D behaves as D\\ cx q 2 , whereas 
its transverse component D± behaves as D± cx q 4 , where 
q is the distance to the boundary.) Our description of 
the chemical reaction problem is based on the solutions 
of Eq. ((2J) in different spatio-temporal domains. Know- 
ing (n), one can establish the temporal behavior of N a ^. 
The brief style of this letter does not allow us to present 
the complete analysis here. Therefore, below we report 
final results, omitting details of the derivation. To clarify 
the results, we also pay special attention to presenting a 
physical picture of the phenomenon. 

We find that the chemical reaction (which starts at 
t = 0) undergoes the following four stages: 

I. Formation of stripes in the bulk. Advection cre- 
ates from an initially smooth distribution a striped struc- 
ture of alternating domains of A and B Q . The stripes 
become dynamically thinner, i.e. inhomogeneities of 
smaller and smaller scales are produced. Once the width 
of the stripe decreases down to the diffusive scale, r*d, 
the stripe collapses (wiped out by the diffusion-limited 
chemical reaction) in a time ~ A . Since the stretching 
(contraction) process leadin g to creation of the stripes is 
exponential in time 0, S H HH EE El , the initial stage 
(when the r^-stripes are formed) lasts for n ~ ln(Sc)/A, 
i.e. just the time required for the cascade of passive scalar 
to run from L down scale to r^. Even though the inter- 
facial area increases exponentially during the first stage, 
N a ,b do not vary significantly. By the end of this stage 
the bulk parts of N a $ begin to decay rapidly (exponen- 
tially), with a decrement of the order of A, according to 
the law of the passive scalar decay in an unbounded spa- 
tially smooth flow EjIEJ- Thus, after the first stage the 
chemicals remain mainly in the peripheral region. 

Notice also, that after the first stage, stripes of differ- 
ent widths, distributed between r<2 and L, are present in 
the bulk. (This multi-scale structure is also seen in the 
passive scalar deca y ex periment EiEa, and in binary- 
reaction numerics When the r^-wide stripe, say, 
of the chemical A collapses, then two nearby stripes of 
the chemical B form one wider stripe. Thus, collapse 
of r^-narrow stripes is accompanied by creation of wider 
stripes, which are shrunk by the flow in turn, and so on 
and so forth. 

II. Peripheral-region- dominated dynamics. The same 
process of layered structure formation takes place in the 
peripheral domain as well. However, advection, which is 
statistically isotropic in the bulk, is strongly anisotropic 



in the peripheral domain, where advection is more effi- 
cient in the direction along the boundary than in the nor- 
mal direction. This anisotropy causes the layers in the 
peripheral domain to stretch mainly along the boundary. 
The stripes closer to the boundary shrink slower than the 
remote ones, since the normal to the boundary compo- 
nent of the stretching rate decreases as one approaches 
the boundary. Therefore the developed layered structure 
(i.e. the one which contains stripes of the diffusive scale 
width) occupies a part of the peripheral region where the 
amounts of A and B become negligible. Thus, the empty 
of chemicals region, formed in the bulk by the end of the 
first stage, starts to expand towards the boundary. As 
a result, the chemicals are arranged mainly within a 5- 
vicinity of the boundary (wall), S ~ L/vXt, where the 
concentrations of the chemicals remain practically un- 
changed. Outside this layer, at L ^> q ^> S (where q 
is the separation from the boundary), the concentration 
of chemicals decreases algebraically (n a j) oc t~ 3 / 2 q~ 3 . 
During this stage the overall amounts of chemicals de- 
crease as S(t), that is cx \j\fi. The spatio-temporal pic- 
ture explained above follows from the universal form of 
the velocity field profile in the proximity of the boundary. 
This stage lasts for r 2 ~ y/Sc/X, i.e. until S shrinks to 
the width of the boundary layer, ru- 

III. Boundary-layer-dominated dynamics. Chemicals 
remain mainly within the r^-thin (not varying with time) 
vicinity of the boundary. The boundary layer width, r&j , 
is still much larger than the reaction zone size (defined for 
the boundary region), so that the passive scalar descrip- 
tion applies. The interfacial area where the chemicals 
interact does not change significantly anymore. Thus, 
due to linear relation between flux of chemicals to the 
interface and their concentrations, the algebraic decay 
switches to an exponential one, i.e. (n a ,b) oc exp(— jt), 
for t ^> r 2 , where 7 ~ \/\/~Sc ~ L _1 -\/Ak. (Note, also, 
that the slow-exponential regime, derived in |l8j for the 
passive scalar, is consistent with the experimental obser- 
vations of 14].) Then, N a .b(t) cx exp(— jt). Chemicals 
are mainly concentrated inside the diffusion boundary 
layer. Outside the boundary layer (at q ^> ru) one of 
the chemicals prevails and its concentration decays alge- 
braically, cx 1/q 3 . The passive scalar description in the 
vicinity of the boundary layer is broken when r c h , which 
grows exponentially with time, becomes of the order of 
ru, i.e. when at the boundary n a ,b becomes ~ A/ (Ry/Sc). 
One concludes that the duration of the boundary-layer- 
dominated stage is r 3 ~ 7" 1 \n(RnoVSc/X), where n is 
the initial concentration of the chemicals. 

IV. Nonlinear stage. By the end of the previous 
stage, advection and diffusion homogenize the remain- 
ing amounts of the chemicals, first, within the bound- 
ary layer and later over the entire reservoir. After that 
there are no inhomogeneities of n a j, left in the sys- 
tem. A purely homogeneous kinetic process takes over: 
diV Q) t,/di = —RVN a Ni, (where V is the chemical reactor 
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volume). Thus, iV 0) & oc 1/t, during the final stage. 

If N a 7^ Nb then the proposed scheme is valid until 
N ai b become of the order of \N a — N b \. Then N a satu- 
rates to a constant (if N a > N b ) and N b disappears expo- 
nentially, oc exp(— RN a t/V). (Note that the exponential 
decay starts after a short intermediate stage character- 
ized by complete homogenization of A due to advection 
and diffusion.) 

Let us now discuss the effect of unequal diffusion co- 
efficients, still assuming that J K a ,b/^ <C L. If Da ^S> 1 
then during the first stages, the chemical length r c h is 
(as above) much smaller than all other scales. This 
problem can also be reduced to a linear one consider- 
ing the advection-diffusion equations in domains popu- 
lated by different species, supplemented by the condition 
that fluxes of the two chemicals towards the interface are 
equal. During the first two stages, the evolution is con- 
trolled by the stripe formation process which is insensi- 
tive to the diffusion. During the latter, third and fourth, 
stages of the evolution in the uneven K a ~ n b case the 
chemicals evolve similarly to what was described above 
for the. Thus, the above description ap- 

plies to the general, K a ^ Kb, case as well. 

We conclude with some general remarks. A compli- 
cated spatio-temporal behavior for the binary chemical 
reaction in a chaotic flow is established. Evolution of 
the chemicals near the boundary (where mixing is slower 
than in the bulk) determines the intermediate stages of 
the reaction. Those boundary-dominated stages were not 
singled out in the previous publications on the subject 
simply because the stages are not observed in an infi- 
nite |j| or periodic [19( flow systems. In our setting the 
lamellar structure (which is statistically isotropic in bulk 
and strongly anisotropic near boundaries) is dynamically 
generated by advection. (This situation is essentially dif- 
ferent from the one considered in j^j, where the lamel- 
lar structure is created initially and no advection par- 
ticipates in subsequent evolution.) We focused here on 
large scale chaotic flows with the size of the box being 
of the order of the major scale of the flow. However, 
it is also of interest for applications to describe chem- 
ical reaction acceleration in turbulent flows, which are 
smooth only inside the viscous range of scales. In this 
case with a large value of the viscous to dissipative scales 
ratio, a consideration similar to those presented in the 
letter is applicable. We plan to examine the more com- 
plicated case in the future. Also, the approach developed 
in the paper is generalizable for other, more complicated, 
types of chemical reaction, e.g. competing chemical re- 
actions. For completeness, let us also mention another 
case of interest which is realized at moderate Da, large 
Sc and if one of chemicals is present in abundance. The 
joint effect of advection and chemistry is different in this 
case (than in the problem discussed in this letter), even 
though rich multi-scale structure of spatial correlations 



is also revealed 21]. A final remark concerns the validity 
of the hydrodynamic description of the chemical reac- 
tion dynamics. It is known that the character of spatial 
fluctuations in the initial distribution of chemicals may 
essentially influence the long-time behavior in diffusion- 
limited chemical systems [23, H3 ■ In some cases (of low 
space dimensionality, d < 2) large scale renormalization 
of the concentration fields due to the small scale fluctu- 
ations could be important. (See, e.g., [13|.) In our case, 
however, this does not happen because the long-time cor- 
relations are completely destroyed by chaotic advection. 

We thank J. Cardy, G. D. Doolen, A. Fouxon, A. Gro- 
isman, I. Kolokolov, V. Steinberg, P. Tabeling, and Z. 
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